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Research Article 

Proteomic signature of reversine-treated 
murine fibroblasts by 2-D difference gel 
electrophoresis and MS: Possible 
associations with cell signalling networks 

Recent advances in stem cell biology have demonstrated that terminally differentiated 
adult cells can be induced to de-differentiate into progenitor cells (induced stem cells) 
upon proper stimuli. This has been achieved by the induced expression of key regulatory 
genes by retro- or lenti-viral systems. On the other hand, synthetic "small molecules" can 
also induce de-differentiation and may represent a potentially safer approach as 
compared with genetic manipulation. Along this line, a synthetic purine called "rever- 
sine" has been shown to induce the de-differentiation of fibroblasts into mesenchymal 
stem -cell-like progenitors, which can be successively induced to differentiate into skeletal 
musde, smooth muscle and bone cells. The mechanism whereby reversine is able to 
achieve de-differentiation has yet to be clarified. In this context, we defined the protein 
changes induced by reversine treatment in murine fibroblasts by 2-D difference gel 
electrophoresis, coupled with MS. Proteins involved in cytoskeletal and cell shape 
remodeling, RNA export, degradation, folding, stress control and ATP production were 
found to be remarkably changed after reversine treatment. Ingenuity pathway analysis 
(I PA) predicted that these protein pattern changes enabled to propose that about 40 
proteins might be associated to several biological functional networks, including cellular 
assembly, cell signaling and cell death. Altogether our data confirm the intrinsic 
complexity of the de-differentiation process induced by reversine and suggest more 
selected approaches to investigate the action mechanism of this small molecule. 
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1 Introduction 

Over the past few years, a new discovery has generated wider 
perspectives in stem cell research. In fact, it was shown that 
terminally differentiated fibroblasts can be re-programmed 
into progenitor cells, almost indistinguishable fi-om bona fide 
embryonic stem cells, as they possess almost identical gene 
expression profiles and esdiibit analogous differentiation 
potential, including teratoma formation capability. In parti- 
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cular, S. Yamanaka's group found that the over-expression of 
four transcription factors (Oct 3/4, Sox2, c-Myc and Klf4) is 
able to reprogram somatic cells, namely skin fibroblasts, into 
embryonic-like stem cells, which were called induced 
pluripotent stem cells [1]. These findings generated great 
interest in the scientific community although the manipula- 
tion of the cell genome with the use of retro- or lenti-viruses 
rose some concerns about the safety of this experimental 
approach. In parallel, small molecules were also shown to be 
able to induce the de-differentiation process [2]. A synthetic 
purine, called "reversine", has been shown to induce the de- 
differentiation of myoblasts into mesenchymal progenitors 
[3], and reversine-treated human or murine fibroblasts could 
be induced to differentiate into skeletal muscle, smooth 
muscle and bone [4]. Thus, it might be envisioned that the 
use of small molecules (such as reversine) to increase cellular 
plasticity (that is to prompt terminally differentiated cells, 
such as fibroblasts, to trans-differentiate into different cell 
types) could be a possibly safer alternative to the genetic 
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manipulation approach, and would also provide an additional 
tool for investigating the complex mechanisms regulating cell 
differentiation. Along this line, an initial study on reversine 
mechanism of action showed that the molecule inhibits non- 
muscle myosin II heavy chain, and MAP kinase and that both 
activities are required for reversine effect [5]. Moreover, it was 
reported that some polycomb genes, which are normally 
involved in chromatin-based gene silencing, are activated by 
reversine [6]. These and other evidence [7-11] support the 
preliminary notion that reversine induces an initial growth 
arrest, which seems to be a key step for the reprogramming 
process. In order to fully consolidate this hypothesis, it will be 
of interest to define the protein pattern involved in growth 
arrest induction. A previous proteomic investigation on the 
eflFects of reversine on C2C12 (an immortalized carcinogenic 
cell line) based on conventional 2-DE on a pH 4-7 gradient, 
indicated a change in 32 protein spots, 22 of them identified 
[6]. The analysis revealed a differential expression of proteins 
involved in cell cycle control, protein biosynthesis and 
metabolism, cytoskeletal rearrangement, proteasome activity 
and electron transport. Based on these, a proteomic analysis 
on primary cultures of dermal fibroblasts, which represent a 
more attractive cell source as compared with C2C12, seems 
very appropriate, also in view of their potential therapeutic 
application. The aim of the present study was to define 
changes in protein expression induced by reversine treatment 
in primary cultures of murine fibroblasts using a more 
advanced analytical technology that is the 2-D difference gel 
electrophoresis (DIGE) coupled with MS. This technology 
allows a precise detection of minor differences of protein 
expression across multiple samples simultaneously with 
statistical confidence by using DeCyder software [12-14]. 
The presence of an internal standard, containing an aliquot of 
each sample under analysis, allows sample normalization. In 
addition, the linearity of the abundance changes, over four 
orders of magnitude, and tiie reproducibility of tiie metiiod 
provide more reliable results tiian conventional approaches 
based on tiie comparison of tiie brightness of gel images 
obtained by conventional staining [15-19]. Proteins involved 
in cytoskeletal and cell shape remodeling, RNA export, 
degradation machinery, folding, stress control and ATP 
production were found to be changed after 4 days of 5 
reversine treatment. Ingenuity pathway analysis (IPA) 
enabled to predict possible associations of tiiese protein 
changes with some key biological networks. 



2 Materials and methods 
2.1 Cell culture and treatment 

Primary mouse dermal fibroblasts were obtained from 
transgenic mice expressing green fluorescent protein [4]. 
They were cultured in DM EM high glucose, supplemented 
with 4mM glutamine, lOOIU/mL penicillin, lOOng/mL 
streptomycin and 10%v/v PBS. Cultures were performed at 
37°C in a humidified incubator with 5% CO2 and 95% air. 
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2.2 Treatment of fibroblasts with reversine 

Reversine was prepared according to the published proce- 
dure [3] and purity (>98%) was checked by HPLC and LC- 
MS analysis. Mouse fibroblasts (8-10 x 10*) were plated in 
100 mm dishes in DMPM supplemented with 10% PBS. 
Pibroblasts were treated with reversine, dissolved in DMSO, 
at concentration of 5 ^M in 10% PBS DMPM, 18-24 h after 
seeding. Control cells were incubated with 0.05% DMSO. 
Treatment was carried out for 4 days without growth 
medium changes. 



2.3 Cell cycle analysis 

Control and treated cells were harvested and fixed in 70% 
ethanol and kept at 4°C before staining. Pixed cells were 
resuspended in 1 ml of a solution containing 5 ng/mL 
propidium iodide (Sigma) in PBS (Gibco Brl), 25 fiL RNAse 
(Sigma) 1 mg/mL, 25 of Nonidet P40 (Sigma) 0.15% in 
water, and stained overnight at 4°C in tiie dark. Cell analysis 
was performed on at least 20 000 events for each sample by 
FACSCalibur System (BD) and DNA prole was analyzed by 
MODFit 3.0 (Verity Software House). 



2.4 Protein extraction 

Cultured reversine-treated and un-treated cells were washed 
twice with PBS and harvested in ice-cold PBS by scraping. 
After centrifugation at 800 x g for lOmin at 4°C the pellet 
was suspended in lysis buffer (7 M urea, 2 M thiourea, 4% 
CHAPS, 30 mM tris and 1 mM PMSP) and 0.1 mM PMSP 
and solubilized by sonication on ice. Proteins were 
selectively precipitated using PlusOne 2-D-Clean up kit 
(GE Healthcare) in order to remove non-protein impurities 
from samples, and re-suspended in lysis buffer. Protein 
extracts were adjusted to pH 8.5 by addition of 1 M NaOH. 
Protein concentration was determined with PlusOne 2-D- 
Quant kit (GE Healtiicare). 



2.5 Protein labelling and experimental design 

Protein extracts (50 |xg) from each set type of three separated 
culture plates, were labelled with 400pmol Cy5 dye 
(CyDyes, GE Healthcare), whereas the internal standard, 
generated by pooling together an aliquot of each reversine- 
treated fibroblasts extracts and each control, was labelled 
with Cy3 dye (CyDyes, GE Healthcare). The minimal 
labelling was performed according to manufacturer's 
recommendations by incubating samples on ice in the dark 
for 30 min. The labelling reaction was quenched with 1 |j,L l- 
lysine 10 mM on ice for 10 min in the dark. The "two dyes" 
protocol was adopted. This experimental design was 
performed by acquiring the images with a two-laser scanner 
(Typhoon 9200, GE Healthcare) and utilizing the combina- 
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tion of Cy3:Cy5, due to their labelling efficienqr and 
reliability compared with other dye combination, as 
described previously [12-15]. 

No preferential protein labelling or background noise 
were observed, but a larger gel number was required 
compared with the three dyes protocol. By utilizing the two 
dyes combination and the manufacturer's dye/protein ratio, 
no dye swap was needed since all samples that undergo the 
statistical analysis were labelled with the same dye (Cy5) and 
normalized against the same internal standard (labelled 
with Cy3). 



2.6 2-DDIGE 

Before lEF, labelled samples were diluted in an equal 
volume of 2 X sample buffer containing 130 mM DTT and 
2%v/v IPG buffer (GE Healthcare). Individual samples 
(40 |xg) were combined with an equal amount of internal 
standard; rehydration buffer (7 M urea, 2 M thiourea, 2% 
CHAPS, 65 mM DTT, 0.5% IPG buffer pH 3.5-9.5 and BBF 
in traces) was added to a final volume of 450 [iL. Samples 
were separated on 24 cm, pH 3-10 non-linear (NL) gradient 
IPGstrips (GE Healthcare), applying the following multi- 
step lEF protocol: 200V (2h), 500V (2h), 1000V(2h), 
2000V (Ih), 3000V (Ih) and 8000 V until a total of 
64 000Vh was reached, using an IPGphor electrophoresis 
unit (GE Healthcare). Protein extracts from each sample 
were run in duplicate. After focusing IPGstrips were 
equilibrated in an SDS-reducing buffer (6 M urea, 2% 
SDS, 20% glycerol, 375 mM Tris-IICl, pM 8.8, 65 mM DTT) 
for 15 min, then alkylated for 8 min in the same buffer 
containing 135 mM iodoacetamide instead of DTT. Second 
dimension was carried out in 20 x 25 cm, 12%T, 2.5%C, 
over the total of monomers, constant concentration poly- 
acrylamide gels at 20°C and 15 mA per gel using the Ettan 
Dalt II system (GE Healthcare). 



2.7 Image acquisition and analysis 

CyDye-labelled gels were visualized using a Typhoon 9200 
laser scanner (GE Healthcare). Excitation and emission 
wavelengths were chosen according to manufacturer's 

recommendations (532 and 633 nm laser beams; 580 and 
670 nm emission filters). Spot detection was performed 
using DeCyder DIA module V. 6.5 (difference in-gel 
analysis, GE Healthcare). Inter-gel matching and statistical 
analysis were performed using a DeCyder BVA module V. 
6.5 (biological variance analysis, GE Healthcare). For each 
gel undergoing the co-detection procedure, the estimated 
number of spots was set to 10 000, filters parameters were 
set as follows: slope >1.2, minimal area cut-off <250 and 
peak height < 14. Statistically significant differences were 
computed by Student's t-test, the significance level was set at 
p<0.01. False discovery rate was applied as multiple test 
correction in order to keep the overall error rate as lower as 



possible. Only proteins with spot volumes consistently 
different in all replicates were considered differentially 
expressed. 



2.8 Protein identification by MALDI-TOF and ESI-MS 

For protein identification, in order to provide a proper 
amount of peptides for MS analysis, semi-preparative gels, 
containing 400 ng of total protein extract per strip, were 
loaded with unlabelled sample; electrophoretic conditions 
were the same as 2-D DICE, except that gels were stained 
with a protein fluorescent stain, as recommended by 
manufacturer (Deep purple, 5 mL/L, GE Healthcare). 
Images acquisition was performed using Typhoon 9200 
laser scanner. Spots of interest were excised from gel using 
the Ettan spot picker robotic system (GE Healthcare), 
destained in 50% methanol/50 mM ammonium bicarbonate 
and incubated with 30 of 4ng/nL trypsin (Promega) 
dissolved in 10 mM ammonium bicarbonate for 16 h at 
37°C. Released peptides were subjected to reverse phase 
chromatography (Zip-Tip C18 micro, Millipore), eluted with 
50% acetonitrile/1% formic acid. An aliquot of 0.35 |xL of 
peptides mixture was spotted onto the sample plate of an 
Ettan MALDI-TOF Pro (GE Healthcare) mass spectrometer; 
an equal volume of lOmg/mL CHCA matrix dissolved in 
70% acetonitrile/30% 50 mM citric acid was applied and 
spots were air dried at room temperature. MS proceeded at 
an accelerating voltage of 20 kV and spectra were externally 
calibrated using Peptide Mix 4 calibration mixture (Laserbio 
Labs); 256 laser shots were taken per spectrum. Proteins 
were identified by comparing the digest peaks with a 
computer-generated database of tryptic peptides from 
known proteins using the built-in search engine Profound. 
Search was carried out by correlation of uninterpreted 
spectra to Rodentia entries in NCBInr database. For 
confirmation of identifications spectra were also searched 
against this database using MASCOT, which utilizes a 
robust probabilistic scoring algorithm. One missed cleavage 
per peptide was allowed and carbamidomethylation, as fixed 
modification, and methionine oxidation, as variable modi- 
fication, were set. Peptide mass tolerance was set at 50ppm. 
In cases where this approach was unsuccessfijl, additional 
searches were performed using ESI-MS/MS. MS/MS were 
recorded using a HCT Ultra mass spectrometer (Bruker 
Daltonics) interfaced to a MDLC capillary chromatograph 
(GE Healthcare). The samples were dissolved in 0.1% 
aqueous formic acid, injected onto a 0.075 x 150 mm Zorbax 
300SB-C18 column (Agilent Technologies), and eluted with 
an acetonitrile/0.1% formic acid gradient. The capillary 
voltage was set to -1600 V, and data-dependent MS/MS 
acquisitions were performed on precursors with charge 
states of 2, 3 or 4 over a survey mass range of 300-1500. The 
collision gas was helium. Proteins were identified by 
correlation of uninterpreted MS/MS to Mus musculus entries 
in NCBInr, using MASCOT software. No mass and pi 
constraints were applied. One missed cleavage per peptide 
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was allowed, and the fragment ion tolerance window was set 
to 0.3 Da. Carbamidomethylation of cysteine was set as fixed 
modification, whereas methionine oxidation as variable 
modification. Where identification was based on a single 
matching peptide, the sequence was confirmed by DeNovo 
Sequencing (Biotools, Bruker Daltonics). 



2.9 SDS Electrophoresis for HMW protein separation 
and quantitation 

SDS electrophoresis was performed on controls and 
reversine-treated cellular extracts on a 37% glycerol, 6% 
polyacrilamide gel according to [20]. Samples (50 \ig) were 
run in duplicate at 100 V, overnight. Gels were stained with 
SYPRO Orange (Molecular Probes) and scanned using a 
570 nm emission filter on Typhoon 9200 laser scanner (GE 
Healthcare). Protein band quantification was achieved using 
Image Quant software (Molecular Dynamics). Molecular 
weight markers containing MHC from rabbit (212 kDa) and 
bovine a2-macroglobulin (170 kDa) were run in a separate 
lane. Differences between groups were computed by a two- 
tail t-test and setting the significativity level at p<0.05. 
Proteins over 200 kDa, differentially expressed or not, were 
excised from gel and hydrolyscd by trypsin, as previously 
described, after reduction and alkylation of cysteine residues 
by 10 mM DTT (100 |xL, 45 min at 56 'C) and 55 mM 
iodoacetamide (100 \iL, 30 min at room temperature in the 
dark) respectively. Proteins were identified by ESI-MS/MS. 



2.10 Immunofluorescence 

Cells were washed with PBS, fixed with 4%w/v parafor- 
maldehyde in PBS and permeabilized with 0.1% w/v Triton 
X-100 in PBS with l%w/v bovine serum albumin for 
30 min. Then, cells were incubated with goat polyclonal anti- 
lamin A/C antibody (Santa Cruz) (1:50 dilution) for 2h at 
room temperature. After incubation cells were washed three 
times in PBS and incubated with the appropriate TRIC- 
conjugated secondary antibody (Sigma) for Ih at room 
temperature. After washing in PBS, cells were analyzed 
under a fluorescence microscope (1x51 Olympus, equipped 
witii an XCite Series 120 fluorescence apparatus, Hamma- 
matsu C8484 digital camera and CellF Imaging Software). 
Images of reversine-treated and un-treated cells were 
acquired in the same experimental conditions (for TRITC, 
acquisition time 200 ms, gain is 1) to avoid possible artifacts. 
Moreover, a surface plot test was done witii ImageJ software 
to compare the fluorescence signal acquired in the two 
images. 



2.11 Immunoblotting 

To correlate MALDI-TOF results with changes in p/s and 
apparent molecular mass provided by 2-D gels, annexins A3 



were confirmed by immunoblotting. For that, a total of 
100 |j,g of control and reversine-treated protein extracts were 
resolved by 2-D PAGE, transferred and blocked onto a PVDF 
membrane (300 mA; 180 min) utilizing a Transblot Cell 
from GE Healthcare. The membranes were blocked over- 
night in TBS (20 mM Tris, 137 mM NaCl, 0.1% Tween, pH 
7.5) containing 2% Blocking Agent powder (GE Healthcare). 
Subsequently the membranes were incubated with mono- 
clonal anti-annexin A3 (1:30 000 as primary antibody; 
Abnova) and HRP conjugate anti-mouse (1:70 000 as 
secondary antibody; GE Healthcare). Proteins were 
visualized by chemiluminescence using ECL Advance kit 
(GE Healthcare). 

To investigate protein ubiquitination, a 2-D immuno- 
blotting was performed with the use of a highly specific 
primary antibody. One hundred microgram of protein 
extract from control and reversine-treated fibroblasts 
respectively was resolved by 2-D PAGE as described 
previously and transferred onto PVDF membranes. 
Membranes were then incubated with antibodies diluted as 
follows: monoclonal anti-ubiquitin (1:20 000 as primary 
antibody. Cell Signaling Technology) and HRP conjugate 
anti-mouse (1:50000 as secondary antibody, GE Health- 
care). Proteins were visualized by chemiluminescence using 
ECL Advance kit (GE Healthcare). 



2.12 Bioinformatic data analysis 

Differentially expressed proteins were evaluated by I PA 
(Ingenuity System, Montain View, CA, USA). IPA is a 
software application that enables to identify the biological 
mechanisms, pathways and functions matching a particular 
dataset of proteins. IPA is based on a database obtained by 
abstracting and interconnecting a large fraction of the 
biomedical literature according to a very strict algorithm. 
This database integrates protein fimctions, cellular localiza- 
tion, small molecules and disease inter-relationships. The 
networks are displayed graphically as nodes, representing 
individual proteins and edges representing the biological 
relation between nodes. Networks are ordered by score and 
optimized including as many differentially expressed 
proteins as possible. A p-score {i.e. — log(p-value)) for each 
possible network is computed. Therefore, networks with 
scores of 2 or higher have at least 99% confidence of not 
being generated by random chance alone. 



3 Results and discussion 

A number of molecules and genes have been already defined as 
targets of reversine. However, a precise mechanism of the 
reversine-mediated de-diflErentiation process has yet to be 
disclosed. Therefore, a proteomic approach could represent an 
essential tool to provide better insights into reversine's mode of 
action. Proteomic analyses were carried out on approximately 
1 X 10'' reversine-treated and control murine fibroblasts, isolated 
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from transgenic mice expressing green fluorescent protein. 
Fibroblasts were grown for 4 days in the presence of 5 |j.M 
reversine in growth medium as described by Anastasia et al. [4]. 
Control cells were treated with the same volume of DMSO, 
which was used to dissolve reversine. Treatment eflFicacy was 
assessed by cell cycle analysis and flow cytofluorimetry. The 
results indicated the formation of a tetraploid cell population 
already after 3 h in treated cells, reaching about 90% at the end 
of the treatment, as shown in Fig. 1. In addition, 20% of treated 
fibroblasts shifted from a diploid/mononuclear (2n) to a 
tetraploid/binudear (2 x 2n) population, confirming previous 
observations [3]. 

To assess changes of the protein pattern due to the 
reversine treatment, proteins from cell extracts were sepa- 
rated by a 2-D DIGE methodology, utilizing a 24 cm, pH 
3-10 NL gradient as first dimension. This method labels 
proteins with fluorescent dyes before 2-DE, enabling accu- 
rate analysis of differences in protein abundance between 
samples. A representative 2-D map of the obtained fibroblast 
proteome is shown in Fig. 2. Overall, approximately 2700 
spots were matched per gel, and 75 were differentially 
expressed (t-test, p<0.01) in treated versus untreated cells. 
The protein load increment, which is necessary for MS 
identification, induces a loss of isoforms resolution and a 
consequent minor number of represented spots [21]. Thus, 
only 66 of the differentially expressed spots were picked up 
and 53 of them identified by MALDI-TOF and ESI-MS/MS. 
Some of the differentially changed spots were very faint and 
close to the limit of 2-DDIGE, consequently, after gel exci- 
sion and digestion, these spots were below the sensitivity of 
our mass instruments and remained undetectable. Identi- 
fied spots are indicated and numbered in Fig. 2. For each 
spot the statistical analysis (average ratio and p-values), the 
identification parameters and protein accession numbers, 
together with the molecular mass and pi, are listed in 
Table 1. Detailed MS data are listed in Table IS as 
Supporting Information. Average ratios for each differen- 
tially expressed protein isoforms are also reported later on in 
the text to facilitate the possible biological interpretation. 



3.1 Semi-quantitation of the high-molecular-weight 
proteins by SDS electrophoresis 

Figure 3 shows qualitative results on high-molecular-weight 
proteins (over 200 kDa), which could not be detected by 
2-DE analysis, by SDS electrophoresis as described in 
Section 2. Differential analysis, conducted by Image Quant 
after Sypro Orange staining, indicated that only two bands 
over 200 kDa were differentially changed. The two bands 
over 200 kDa were excised, digested with trypsin and 
analyzed by ESI-MS/MS. ESI-MS/MS results are summar- 
ized in Table 2 and in Supporting Information Table IS. 
Two ubiquitously expressed isoforms of Filamin, A and B, 
were significantly (p<0.05) underexpressed in reversine- 
treated samples. Filamins are known to bind actin, different 
molecules and ions involved in the cellular signaling and 
transcriptional regulation [22-24]. In particular Meng et al. 
[25] suggested that filamin A plays a role in mechanisms 
making cells more vulnerable to several DNA-damaging 
agents and that lack of filamin A causes a delay in G2 
recovery after irradiation. 

Myosin II heavy chain isoforms were also identified. No 
changes in protein abundance were visualized (p = 0.09), 
confirming that reversine treatment does not induce an 
increment of non-contractile high-molecular-weight 
myosins II but its role is related to its ATPase activity 
inhibition [5]. To better clarify this issue further investiga- 
tions will be required. 

The results provided by this differential analysis support 
data from previous studies, confirming the effects of rever- 
sine on cell cycle, cell proliferation and cell morphology 
[4-6]. According to gene ontology the differentially expres- 
sed proteins were clustered into four major categories: 
cytoskeletal and cell shape, degradation machinery, stress 
control and metabolic proteins. 

To gain better insights into biological significance of the 
quantitative differential results, differentially expressed 
proteins were analyzed by I PA Software (www.ingenuity.- 
com). I PA is a software application that enables to identify 
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Figure 1. Effects of re 



blasts cell cycle. Flow 
cytometry analysis on five 
samples performed after 
72 h of DMSO (control) or 

5nM n 
After n 
for 72 h, the percentage of 
diploid population was 
12.76 -F/-1.1, whereas the 
percentage of tetraploid 
population was 87.24 
+/-5.2. In control samples 
the percentage of diploid/ 
tetraploid population was 
99 +/-^ and 1 -t-/-1, 
respectively. 
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Figure 2. Representative 2-D DIGE map of fibroblast protein 
extract labelled with Cy3. Fifty micrograms of protein extract 
were separated in a NL pH 3-10 IPG strip in the first dimension 
and SDS gel {12% T, 2.5% C) in the second and acquired with a 
532 nm laser beam and 580 nm emission filter. DeCyder software 
allowed to detect approximately 2700 spots. Differential analysis 
of reversine-treated fibroblast cells compared to control {statis- 
tical analysis was performed by BVA module of DeCyder) 
revealed 75 differentially expressed {p<0.01) protein spots. 
Identified proteins were indicated by circles and numbered as 
in Table 1. DeCyder Software 3-D views for differentially 
changed spots 2, 3, 4, 8 and 9 are also showed. 



the biological mechanisms, pathways and functions 
matching a particular dataset of proteins. In this study 55 
differently abundant proteins were uploaded utilizing IPA 
bioinformatic tool. Among them, 42 proteins were involved 
in networks on the basis of IPA criteria. As summarized in 
Table 3, IPA indicated a total of five networks with the 
following scores: 33 (16 molecules involved), 23 (12 mole- 
cules involved), 22 (12 molecules involved), 2 (1 molecule 
involved) and 2 (1 molecule involved), respectively. Two out 
of five networks (1 and 2) shown in Fig. 4A and B, respec- 
tively, appear to overlap. The networks with high score are 
associated with cellular assembly and organization, cancer, 
cell-to-cell signaling and interaction, cell death, neurological 
disease, molecular transport and cellular function and 
maintenance. In particular, proteins were grouped into 
functional categories based on their biochemical role. 



3.2 Cytoskeletal and cell shape remodeling 

Active remodeling of the cytoskeleton is fundamental to 
most biological processes that involve motility and changes 
in cell morphology, including cell migration, spreading, 
division and cell fusion. Changes in cell morphology can be 
elucidated by observing variations of proteins involved in 
actin filaments dynamics and remodeling [26]. Differential 
analysis showed a positive increment of cofilin 1 and Rho 
GDI a ( + 1.3- and + 1.49-fold, respectively) in reversine- 



treated cells. Cofilin 1 promotes actin filament disassembly 
[27], which is present in myotubes and its phospho/ 
dcphosphorylation is essential for cytoskeletal rearrange- 
ment [28]. On the other hand, Rho GDI is an inhibitory 
regulator and promotes the reorganization of actin filaments 
with a consequent cytoskeletal remodeling [29]. Tropomyo- 
sin beta, which stabilizes actin in non-muscle cells, and 
Tropomodulin 3, which blocks elongation and depolymer- 
ization of actin filament, were also found increased (-1-1.54- 
and + 1.67-fold, respectively). 

Another strongly increased protein target of reversine 
action was P-tubulin ( + 1.66-fold), a cytoskeletal protein, 
which is highly abundant in microtubules [30] and that plays 
a crucial role in cytoskeletal remodeling. In particular, an 
increment of P-tubulin promotes the assembly and organi- 
zation of the cytoskeleton into 3-D structures by binding the 
vmculin, moesin and the disassembly of actin depolymeriz- 
ing factor/cofilin [31]. In treated protein extracts, vinculin 
(-1-1.33), and two subunits of moesin (-1-1.79- and -1-2.21-fold) 
were more abundant. 

Rho GDI, cofilin and tropomyosin are known regulators 
of actin filament dynamics and membrane extensions. The 
increased levels of Rho GDI-a (+1.49) rely on the dynamic 
reorganization of the actin cytoskeleton, contributing to 
control membrane protrusions associated with apoptosis 
[32] and cell injury [33]. In addition, cofilin is essential for a 
rapid turnover of actin filaments, actin polymerization and 
membrane protrusions, and has the capacity to determine 
the direction of cell motility [34, 35]. By combining these 
results, we can speculate on a positive action of reversine in 
cytoskeletal reorganization, and on its negative modulation 
of the pro-apoptotic formation of protrusions on the cell 
surface. 

In contrast to the increment of some cytoskeletal 
proteins, nucleophosmin (NPM) was decreased (—1.85-fold). 
This protein is also involved in cell-shape changes, as it is a 
highly phosphorylated protein associated with the nucleolar 
ribonucleoprotein structures, which were previously found 
decreased in reversine-treated cells [36]. NPM is involved in 
many cellular processes but its function remains unknown. 
Amin et al. [37] demonstrated that in HeLa cells, NPM 
causes distortion of nucleolar structure leading to a multiple 
micronuclei formation induced by a cytoskeletal rearrange- 
ment. Moreover a decrement of NPM produces reduced 
cellular proliferation and apoptosis [38]. 

A similar decrement (—1.78-fold) was also observed in a 
tumor protein, translationally controlled 1 (TPTC 1), an anti- 
apoptotic protein involved in calcium binding and micro- 
tubule stabilization [39]. 

Interestingly lamin C, which determines nuclear shape 
and size, was increased ( + 2.41-fold). 

Physiologically the presence of lamin C is related to 
cellular differentiation [40]. It is known that nuclear lamin 
A/C deficiency induces defects in cell mechanics, polariza- 
tion and migration. Actually, lamin A/C deficiency leads to 
nuclear fragility and cytoplasmic mechanical rigidity by 
reducing its elasticity and viscosity [41]. Moreover, lamin 
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cathepsin D. 

It is known that in eucariotic cells the degradation 
machinery can follow lysosome and/or proteasome path- 
ways: both of them were shown here to be deregulated by 
reversine treatment. Proteasome 28 subunit a is a 249 
amino acids small activator of a peptidase, activating 
hydrolysis of small non-ubiquitinated peptides by means of 
the 20S proteasome [43]. However, even if its precise cellular 
fiinction remains only partially imderstood [44] we found its 
expression increased (+1.25). On the contrary, the Protea- 
some subunit a type 6 expression was slightly decreased 
(—1.13). This protein is one of the subunits constituting the 
20S proteasome core, which is responsible for the hydrolitic 
activity. In particular, the a subunit represents one of the 
two outer rings whose function is to bind the 19S regulatory 
subunit, forming a channel essential for the degradation of 
denatured proteins [45]. The effects of proteasome inhibi- 
tion on the stability of various cell cycle regulatory proteins 
are well established and include the arrest of mitogen- 
induced proliferation, inhibition of cell cycle entry and 
apoptosis [46, 47]. 

Together with slight changes on proteasome subunit 
assembly, a differential expression of two enzymes involved 

El ( + 1.74-fold in spot 3 and + 1.79-fold in spot 2) and 
ubiquitin-conjugating enzyme E2 variant 1 (—1.46-fold) was 
observed. The first enzyme plays a role in recruiting 
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Name p/ Molecular weight (Da) Swiss-Prot MASCOT Score Identified peptides Coverage (%) 

accession number 



Filamin-A 5.68 281 193.50 Q8BTM8 1789 71 38 

Filamin-B 5.44 277 752.59 Q80X90 2133 70 39 

Myosin, heavy polypeptide 9, 5.54 226372.08 Q8VDD5 829 25 19 
non-muscle isoform 1 



Table 3. The high-scored biological networks in murine reversine-treated fibroblast cells^' 
Network Genes in network Top functions 



1 ACADL, AC0T2, ALDH2, ATF4, ATP6V1B2, CDK2, CREB1, CTSD, EGF, Cell-to-cell signaling and interaction, endocrine 33 16 
EX0SC6,GSTM1 (includes EG:2944),HNF4 alpha dimer,HSP90B1,HTT, system development and function, molecular 

IL4, IL6, IL12, KITLG (includes EG:4254), LY6E, NFE2L2, NPM1 (includes transport 
EG:18148), PCK2, PSMA6, PSMB7, PSMB8, PSMEl, PSME2, PYCRl, 
retinolc acid, RNF128, SHIVrT2, SWAP70, TPT1 (includes EG22070), 
TUBB 

2 ANXA5, ARHGDIA, ARHGEF1, C190RF10, CFL1, DDAH2, Erm, EZR, F Cell morphology, cellular assembly and organi- 23 12 
Actin, F2RL1, FOS, Ige, LIMK2, LMNA, LY6E, MARCKSLl, MKLl, lUSN, zation, cellular function and maintenance 

IVIYH7, PFN, PIP5K1B, RDX, retinoic acid, RHOA, Rock,TGFBl, TM0D3, 
TNNT2, TPI1, TPM2, Tropomyosin, TTN, UBA1, VASP, VCL 

3 AK1, ALDGA, ANXA3, ANXA7, BAG3, Ca2+, CCT2, CCT3, CCT4, CCT5, Cancer, cell death, neurological Disease 22 12 
CCT7, CCT8, CCT6A, C0L3A1, DNAJ, DNAJC13, ERP29, GLUDl, GSTMl 

(includes EG:14862), HSPA8, HSPA9, IGF1R, INS1, KCNJ11, KLP4, l- 
glutamic acid, SERPINHl, SRI, ST13, STIPl, TCPl, TOPI, TP53, 
UBE2V1, VHL 

4 ETFA, HEXB Carbohydrate metabolism, lipid metabolism, 2 1 

small molecule biochemistry 

5 Selenide, water dikinase, SEPHS1 - 2 1 



a) The focus proteins are indicated with gene names and shown in bold. A score of 2 was considered significant (p<0.01). 



ubiquitin, which is transferred to a target protein through 
the action of the E2 variant 1 enzyme. This step of the 
proteasome machinery seems to be inhibited by reversine 
treatment. 

Beside the degradation via proteasome, a differential 
expression of proteins involved in lysosomal degradation 
was also observed. This study showed a decrement of 
vacuolar adenosinetriphosphatase subunit B (—1.68-fold) 
and an increment of Cathepsin D (-1-1.79-fold). The latter is 
a lysosomal aspartic endoprotease, which plays a crucial role 
in the regulation of apoptosis [48]. 



3.4 Protein folding 

2-D DIGE showed a deregulation (—1.23-fold) of chaperonin 
containing TCP-1 6 subunit involved in folding, particularly 
of actin and tubulin both in vivo and in vitro. Its main 
function is to bind denaturated/neo-synthetized proteins 
and to mantain cellular homeostasis [49, 50]. Coghlin et al. 
[51] pointed out that TCP-1 is not only involved in cell 
growth but is also highly over-expressed in the S -phase of 
cell cycle. 



Additional changes were observed in HSPs whose major 
role is to assist protein folding. In this study we found that 
three different HSPs were differentially expressed: HSP9A, 
47kDa HSP (-1-1.23- and -1-1.66-fold-change, respectively) 
and heat shock cognate 71 kDa protein (HSPaS; -1.47-fold 
change). HSP9A and HSPa8 are members of HSP70 family 
[52] involved in many biological processes such as stress 
response, cell death/proliferation and trafficking control 
[53]. In contrast, tumour rejection antigen gp96, a glyco- 
protein belonging to HSPs' family, that plays a role in 
different cellular functions, such as homeostasis and cell 
proliferation [54], was found decreased (—1.45-fold). More- 
over we found an increment of endoplasmic reticulum 
protein ERp29 precursor (-1-1. 31 -fold). This protein partici- 
pates in protein folding in endoplasmatic reticulum and it 
plays an important role in the processing of secretory 
proteins [55]. 



3.5 Stress control 

A strong increment of two cytosolic proteins involved in 
oxidative stress such as glutathione S-transferase, mu 1 
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4 




Figure 4. Signaling networks 1 (A) and 2 (B) of murine fibroblast 
proteins regulated by reversine treatment. Identified differen- 
tially changed proteins (indicated by their gene names) were 
uploaded in IPA. Edges containing single line indicate protein 
associaiions, arrows indicate expression-conirollmg proiems, 
whereas grey nodes represent the Identified proteins (focus 
proteins). Nodes shapes (C). 



( + 2.07-fold) and mu 2 ( + 2.69-fold) subunits was observed. 
This increment reflects the response to intracellular 
oxidative stress induced by reversine. Dandrea et al. [56] 
suggested that transient S-glutathionylation of protein is a 
cross-protection against morphological changes generated by 
the induction of an oxidative DNA damage and by the 



presence of agents inducing cell toxicity. Moreover, the most 
important role of the glutathione S-transferase enzymes is 
the inactivation of xenobiotic agents [57]. Particularly the mu 
class is active in detoxification of electrophilic compounds, 
including carcinogens, therapeutic drugs, environmental 
toxins and products of oxidative stress, by conjugating 
them with glutathione. Pyrroline-5-carboxylate reductase 1 
(+3.96-fold), a proline biosynthetic enzyme, was also up- 
regulated. Accumulation of proline in cytoplasm is an 
adaptative response to stress. Krishnan et al, showed that 
proline influences glutathione levels during H2O2 stress 
treatment [58]. The level of peroxiredoxin 3 was also changed 
(-1-1.41-fold). Peroxiredoxins are small proteins often identi- 
fied in proteomic studies [59], involved in processes such as 
proliferation, transformation and apoptosis. 

Dimethylarginine dimethylaminohydrolase 2 which 
hydrolyzes N(G),N(G)-dimethyl-L-arginine and N(G)-mono- 
methyl-L-arginine acting as a NOS inhibitor and playing a 
role in nitric oxide generation was also slightly increased 
(-1-1.15-fold), whereas the selenophosphate synthetase 1 was 
slightly decreased (—1.1 3 -fold). This protein is involved in 
ATP binding and selenocysteine synthesis [60]. In addition, 
glutamate dehydrogenase 1 mitochondrial precursor, which 
is a reducing enzyme involved in energy storage, was found 
increased (1.41-fold). 



3.6 Metabolism 

Phosphoenol pymvate carboxy kinase 2, involved in 
gluconeogenesis and nucleotide binding, and fructose- 
bisphosphate aldolase catalyzing the reversible conversion 
of Fru-l,6-biphosphate to glyceraldehyde 3-phosphate and 
dihydroxy-acetone were more abundant ( + 2.03- and +3. 55- 
fold, respectively) in reversine treated cells. Both enzymes, 
however, are gluconeogenic intermediates and elevate the 
intracellular levels of ATP [61, 62] supporting energy 
demand. In contrast the conversion of glyceraldehyde 
3-phosphate to dihydroxy-acetone phosphate was reduced 
due to a decrement of triosephosphate isomerase (-1.54- 
fold). Adenylate kinase 1, a small ubiquitous enzyme that is 
essential for maintenance and cell growth, and the 
mitochondrial serine hydroxymethyltransferase 2, which 
regulates the inter-conversion of serine and glycine to 
generate pyruvate, were increased (-1-1.58- and -1-1.63-fold, 
respectively). 

Five enzymes related to mitochondrial activity were also 
more abundant; Acyl-CoA thioesterase 2, the long chain- 
specific Acyl CoA dehydrogenase, the electron transfer 
flavoprotein subunit a, enoyl CoA hydratase and NADH 
dehydrogenase (ubiquinone) 1 a subcomplex 10 ( + 2.35, + 
1.69, +2.7, +1.53 and +1.63, respectively), suggesting an 
activation of the aerobic metabolism. The acyl-CoA thioes- 
terases are a group of enzymes that catalyze the hydrolysis 
of acyl-CoA to free fatty acid and coenzyme A (CoASH) by 
interacting with long chain-specific acyl-CoA dehy- 
drogenase, providing the potential to regulate intracellular 
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levels of acyl-CoA, free fatty acids and CoASH. The third 
enzyme acts as a specific electron acceptor of dehy- 
drogenases, including five acyl-CoA dehydrogenases trans- 
ferring electrons to the mam mitochondrial respiratory 
chain via electron transfer to ubiquinone oxidoreductase. 
Enoyl-CoA hydratase is involved m the second step of 
P-oxidation of fatty acids by catalyzing the reversible 



confirmcxl by our analysis (+2.1 4-fold). The; mammalian cyto- 
chrome P450 IS a component of membrane-bound mono- 
oxygenase family involved m xenobiotic agents, drug 
metabolism and signalling of endogenous molecules [65]. This 
family of enzymes can be regulated by a number of post- 
translational events m vivo, including availability of catalytic 
treatment, post-translational m<Kiifications and protein-protein 




pH 10 



Figure 5. 2-D DICE of control sample (Cy5 labelled) (A) and close-up of the 45l<Da region containing annexin A3 (spots 27 and 30 in 
Table 1) (B), annexin A3 Identification by monoclonal antibody after electrotransfer to PVDF membrane (C), annexin A3 localization 
In reverslne-treated sample (D) and after antibody recognition (E). Proteins were separated In a NL pH 3-10 IPG strip In the first 
dimension and SDS gel {12% T, 2.5% C) In the second and electrotransferred onto PVDF membranes. The membranes were Incubated 
with monoclonal antl-annexin A3 (1:30000 as primary antibody) and HRP conjugate anti-mouse (1:70000 as secondary antibody). 
Proteins were visualized by chemilumlnescence using ECL Advance kit. 



hydration of trans-2-enoyl-CoA thioesters to the corre- 
sponding 3-hydroxyacyl-CoA thioesters [63], whereas NADH 
dehydrogenase (ubiquinone) 1 a subcomplex 10, is a key 
enzyme of the respiratory chain complex 1. 




3.7 Others 

Among the more abundant proteins we can indicate 
annexins (A3, AS and A7) and calreticulin. The latter is a 
multifunctional protein involved, as annexin, in Ca^"^ 
binding and folding [64] and was -1-1.17-fold increased. 

Annexins 3 and 5 were identified with diflFerent isoforms in 
a 2-D map, whereas annexin 7 was present as a single increased 
spot (-1-1.63). These cytosoUc proteins are involved in binding of 
Ca^"*" enabling the interaction with membrane phospholipids, 
and in membrane-binding regulating cell signaling and prolif- 
eration. Both isoforms of annexins AS were more abundant in 
reversine-treated fibroblasts (+1.13- and +1.46), whereas 
s A3 showed a differential trend of the two isoforms: 
1 3 corresponding to spot 30 was decreased (—1.36), 
whereas isoform of spot 27 was increased (+1.39). Furthermore, 
the representative 2-D map indicates both annexins A3 have a 
more alkaline isoelectric point in comparison with the theore- 
tical one suggesting the possible presence of a post-translational 
modification, the investigation of which is out of the topic. 

Previous proteomic investigations on myoblasts [6], 
indicating an increment of cytochrome P4S0 after 1- 6 n 
treatment (family 19, subfamily a, polypeptide 1) ' 




e 6. Close-up of a 2-D gel for the assessment of protein 
ubiqultlnatlon after Immunoblotting of 100|j.g protein extract 
from control and reverslne-treated fibroblasts. Proteins were 
separated In a NL pH 3-10 IPG strip In the first dimension and 
SDS gel (12% T, 2.5% C) In the second and electrotransferred 
onto PVDF membranes. The membranes were incubated with 
antibodies diluted in the following way: monoclonal anti- 
ubiquitin (1:20000 as primary antibody) and HRP conjugate 
anti-mouse (1:50 000 as secondary antibody). Proteins were 
visualized by chemiluminescence using the ECL Advance l<it. 
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interactions [66, 67]. Particularly, P450 19A1 has been related to 
breast cancer due to its involvement in estrogens biosynthesis. 
The enhanced expression of P450 19A1 actually increases 
estrogen production, resulting in a positive feedback loop to 
drive hormone-responsive breast cancer cell growth [68]. In our 
cellular model we can speculate on the positive role of this over- 
expression. Estrogens and their receptors, play a role in regu- 
lating processes related to start, progression and exit from the 
Gl phase of tiie cell cyde [69, 70]. 

Mtr3 (mRNA transport regulator 3)-homolog, a 
component of the exosome 3'-> 5' exoribonuclease complex 
degrading unstable mRNAs and containing AU-rich 
elements required for the 3'-processing of the 7S pre-RNA to 
the mature 5.8S rRNA, was increased (+ 1.64-fold). 

Another protein involved in m-RNA processing was 
found deregulated: the heterogeneous nuclear ribonucleo- 
protein A2/B1 isoform 2 (—1.16). This protein is involved in 
pre-mRNA processing and its change is a signal (together 
with ubiquitin-conjugating enzyme E2 decrease) of DNA 
damage and transcription detuning [71]. DEAD (Asp-Glu- 
Ala-Asp) box polypeptide 39, involved in pre-mRNA 
splicing and required for the export of mRNAs out 
of the nucleus, was decreased (-1.45-fold). The deregulation 
of these two proteins suggests a direct modulation of 
reversine both in the elimination of the elements, which 
could induce instability, and in the control of the nuclear 
mRNA export. 



Finally we found a variation concerning interleukin-25 
( + 1.41-fold), which is involved in immunity response by 
modulating cytokines expression, such as IL-5 and IL-13; 
however, its regulatory fiinction remains unknown [72]. 



3.8 Validation 

With regard to differential expression of annexin 3 isoforms and 
their localization on the 2-D map, blotting was conducted after 
separation on 24 cm pH 3-10 NL gradient. Figure 5 shows the 
results of detection of the two isoforms in a control and 
reversine-treated samples. Spots are differentially abundant in 
controls versus reversine-treated fibroblasts according to DeCyder 
quantitation trend. In addition, both isoforms were recognized 
suggesting that changes in apparent molecular weight and pJ 
could be a result of a post-translational modification, which will 
be characterized in fijrther studies. 

As to annexin 3, ubiquitinated proteins were also analysed 
by immunoblotting after 2-DE separation under the same 
gradient and electrophoretic conditions. The number and the 
amount of ubiquitinated proteins were lower after reversine 
treatment, supporting the 2-D-DIGE information. Figure 6 
shows a close-up of ubiquitinated proteins id(;ntifi(;d by 
immunoblotting on control and treated sample. 

A fiirther validation of our differential results was 
obtained by quantifying lamin A/C expression and locali- 




murine fibroblasts with 
antibody against lamin A/C 
(red) performed after 72 h 
of DMSO (left panels) or 
5nM reversine treatment 
(right panels). Fluorescence 
microscopy (20 x ) reveals 
higher levels of lamin A/C 
staining In reversine-treated 
cells, as shown In the 
surface plot analysis (lower 
panels). Fluorescence Inten- 
sity scale Is also reported In 
the Inset graphs. 
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zation by immunofluorescence as shown in Fig. 7. The 
results are in agreement with the proteomic differential 
analysis. 



4 Concluding remarks 

Induced stem cells could be a valid source in addition to 
adult and embryonic stem cells in order to obtain the widest 
extensions, in terms of different organs and tissues to be 
replaced in regenerative medicine. The use of "small 
molecules" to produce induced stem cells is being adopted 
by an increasing number of research groups. Therefore, it 
has become crucial to understand the mechanisms under- 
lying the de-differentiation process. This would allow 
finding the ways to increase the efficiency of the process, 
and to assess its safety. In particular, a wide-range 
proteomic investigation coupled with advanced bioinfor- 
matic tools seems a powerful approach to study the effects 
induced by de-differentiating drugs such as reversine. 
Although the raw data output from a comprehensive 
proteomic analysis may result quite difficult to be inter- 
preted, the use of bioinformatic algorithms, as the I PA 
software, may give important hints to recognize or predict 
the major cell pathways affected by reversine treatment. In 
fact, as shown in this study, reversine appears to exert a 
crucial role on specific proteins involved in cell cycle arrest, 
cytoskeletal remodelling, stress response, mRNA export, cell 
signalling and proliferation. These protein changes, as 
anticipated, will require further studies to be validated, but 
constitute an important step-forward towards the under- 
standing of reversine mechanism of action. 
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